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Revivals of the coherent states of a deformed, adiabatically and cyclically varying 
CN ■ oscillator Hamiltonian are examined. The revival time distribution is exactly that 

of Poincare recurrences for a rotation map: only three distinct revival times can 
^\i I occur, with specified weights. A link is thus established between quantum revivals 

and recurrences in a coarse-grained discrete-time dynamical system. 
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Revivals of a non-stationary state of a Hamiltonian may be regarded as the quantum 
analogs of recurrences in a classical dynamical system. In general, when a quantum 
mechanical system is in a state I'?/'), the autocorrelation function \{ilj{0)\ilj{t))\'^ decays from 
its initial value of unity with increasing time t. Under very special circumstances, however, 
it may return to its initial value at certain instants of time, signalling a full revival. 
Typically, however, any return would be to a value in a small neighborhood of unity 
IHI], which we shall term a "near revival". This would be the counterpart of a Poincare 
recurrence in a classical system, although, of course, there is no direct connection 
between the revivals of quantum states and recurrences of the classical phase trajectories 
of the same physical system. Taking into account finite experimental resolution, it is 
evident that near revivals rather than exact revivals are the readily identifiable phenomena 
of interest. Revivals and fractional revivals of wavepackets are of considerable current 
interest in many contexts: atom optics |^ , propagation of coherent light through a Kerr 
medium [Q, Rydberg states 0, interacting atom-radiation (Jaynes-Cummings) systems 
|]1|, H, 1^ with single and multi-photon coherent and squeezed coherent states, quantum 
dynamics of boson systems and spin systems ||^, and so on. 

In this paper, we show that one can indeed establish a close relationship between the 
near revivals of a wavepacket and recurrences in a discrete-time dynamical system, by 
essentially mapping the former problem onto the latter. This leads to a convenient and 
systematic way of analyzing the statistics of revivals. In order to discretize time in an 
unambiguous and natural way, the parameters in the Hamiltonian are varied adiabatically 
and cyclically, with a time period T. Consequently, it is only at instants separated by 
an interval T that the Hamiltonian returns to its original self. It is therefore meaningful 
to look for revivals of a state only at the instants T, 2T, ■ ■ -. In turn, this permits us to 
analyze the problem in terms of recurrences in a discrete-time map. 

As a by-product of the cyclic variation of the parameters in the Hamiltonian, the 
state vector may pick up an extra non-integrable ("geometric") phase [^. One may 
expect the revival times to be affected as a consequence - a priori, to get shifted by the 
time required to cover the angular excess (or decrement) representing the anholonomy in 
the semiclassical limit, namely, the corresponding Hannay angle [l^, 0. This is indeed 



the case [|12[ , and the specific system we work with yields an explicit verification of this 
result as well. For definiteness, we consider a unitarily deformed oscillator Hamiltonian 
whose eigenstates are squeezed generalized coherent states, while the spectrum remains 



linear (equi-spaced). Coherent states and their generalizations |TB| are convenient in this 
regard, as they approximate classical states of radiation. They also provide a setting in 
which anholonomies (non-integrable phases) can be measured experimentally [ll4| , pTSJI and 



analyzed theoretically |TB, T7], in a physical range extending from the extreme quantum 
regime to the semiclassical limit. Although revivals, per se, are of greatest physical interest 
in the case of spectra that are nonlinear (in the quantum numbers), here we restrict our 
attention to a linear spectrum. (In this case exact revivals correspond, as we shall see, 
to a trivial periodicity.) Complications such as fractional revivals |T^ do not arise here. 



and this helps us establish in a clear manner a very interesting link between near revivals, 
anholonomies and recurrences - manifestations, respectively, of quantum interference, 
non-trivial topology in parameter space, and ergodicity in a classical dynamical system. 
We recall that the coherent state \() (C G C) of the standard oscillator Hamiltonian 
H = hu!{a^a+l/2) is given by \() = D{() |0) where |0) is the ground state of H, and D{Q 
is the displacement operator exp {(a^ — (*a). For each value of (, \() is represented in the 
position basis by a Gaussian wavepacket which does not spread with time (under evolution 
governed by H). The centre of the wavepacket oscillates according to the classical laws 
of motion. Using the expansion of |C(0)) = |C) in terms of the eigenstates {|^)} of H, its 
correlation function is easily found to be 

|(C(0)|C(t))r = exp[2|Cp(cosu;t-l)]. (1) 

Thus, in this linear case, revivals of the initial wavepacket simply amount to periodicity 
with period 2tt/uj. To induce a geometric phase, however, we need a Hamiltonian that has 
more than just a single parameter (analogous to u in H) that can be varied. Further, in 
order to have a no n- vanishing anholonomy (Hannay angle) in the semiclassical limit, the 
geometric phase must have a dependence on the quantum number n labelling the states 
|T0| . Both these objectives are achieved if we begin with a deformation of H produced 
by squeezing. For the sake of generality, we include a possible displacement as well, and 
define the transformed Hamiltonian 

H{a,(3) = S{p) D{a) H D\a) S\p) (2) 

where a, /? G C, and S{I3) is the squeezing operator [|19| S{I3) = exp {{jSa^ — (3*0?) /2). 
As D{q) and S{(3) are unitary, H and H are unitarily equivalent and isospectral. H has 
eigenstates |n, a,/?) = S{I3) D{a) \n) { n = 0,1,---), connected by ladder operators a 
and a\ where a = S D a D"^ S^ and [a, a^ = 1. 

Under a cyclic, adiabatic variation of the displacement and squeezing parameters a = 
«! + ia2 and P = Pi + i P2, the state \n , a , P) acquires a Berry phase 7^ that is a 
sum of contributions from the displacement and squeezing parameters, respectively. The 
contribution from the variation of the displacement parameter a can be found from the 
group multiplication law for the elements {D{a)} of the Heisenberg-Weyl group, namely, 

D{a)D{a') = D{a + a') exp [ix{a,a')] (3) 



where x(q^) Q^O is twice the area of the triangle with vertices at 0, a' and a + a' in the 
complex plane. This contribution is independent of n essentially because [a, a^] is just the 
identity operator. On the other hand, the contribution from the variation of the squeezing 
parameter (3 is more involved. The squeezing operator S{I3) represents an element of the 
group 5f/(l,l), with generators K+ = a^^/2, K_ = o? jl and K^ = [K_,K+]/2 = 
{aa'' + a^a)/4. The multiplication rule is of the form [|1^ 

S{(3) S{(3') = S0') exp [i 0(/3, /?') Ko] (4) 

where /? is a certain function of /? and P'. The phase 0(/5, P') is related to the area of a 
certain geodesic triangle on the upper sheet of the hyperboloid x ■ x = Xq — xf — X2 = 1, 
where (setting arg/5 = (p) x = (cosh2|/3|, — sinh2|/3| cos(p, sinh2|/3| simp), on which the 
parameters of SU{1, 1) live. The geometric phase is essentially the solid angle subtended 
at the origin x = by the invariant area enclosed on the hyperboloid [^]. Since x-x = 1, 
this is just the surface integral over {dxi dx2)/xo on the hyperboloid. In terms of P this 
integral is 



B 



/.^,!^. (5) 



the integration running over the area enclosed by the loop traversed in the /5-plane in 
a cyclic variation of the squeezing parameters. Further, owing to the presence of the 
diagonal generator Kq in Eq. (4), as opposed to the unit operator in the case of D{a), 
there is an additional factor {n + |) in the geometric phase. Collecting these results, we 
arrive at the expression 

7, = -2A-{n + ^)B, (6) 

where A is the area of the loop traversed in the a-plane in one cyclic variation of ai and 
tt2. The linear dependence of 7„ on n [|l^ implies that A and B are determined by 70 and 
7i (and vice versa): A = (71 — 870) /4, i3 = 70 — 71. 

We now construct the coherent state \z) in the tilde-basis, i.e., d\z) = z\z) {z G C), 
so that 

/ \z\'^\ °° z"^ 

\z) = exp — tt) Yl —n 1^ ' « ' /^) • C^) 



n=0 Vn! 

In the present instance, an equivalent way of obtaining \z) is by operating on the ground 
state |0, a,/5) by the displacement operator exp(2;a''" — z*a). Consider time evolution 
governed by the Hamiltonian H, while H itself changes adiabatically owing to a cyclic 
variation of the parameters a and /3 with a time period T >> uj^^ . An initial state |;i;(0)) 



\z{T)) = exp I — —j Y, -^ exp ^-i (^7„ - (^ + 2 ) ^^j 



|n,a,/?) . (8) 



= \z) (given by Eq. (7)) evolves at time T to 

|2\ 00 ^n 
. n=0 V^ 

Using the results quoted in Eqs. (5) and (6) for 7^ , the correlation function now works 
out to 

\{z{0)\z{T))\^ = exp[2|z|2(cos(cuT + S)-l)] . (9) 

Therefore revivals occur at times T, 2T, ■ ■ ■ provided uoT + B = 27rp, where p = 0, 1, • • • - 
that is, provided the cyclic variation of the parameters in the Hamiltonian is carried out in 
a time period equal to one of the values (27r—i3)/u;, {4:Tt—B)/uj, ■■■. This is to be compared 
with the original revival times 2iTp/uj (cf. Eq. (1)). As —d'jn/dn = i3 (= 70 — 71 in 
the present instance) is just the Hannay angle [|TU], we have here an explicit verification 



of the staggering of the revival times by precisely the time required to cover this angular 
excess (or decrement, depending on its sign) [jl2|. 



Next, we note that there is a convenient and natural way to describe the evolution of 
the initial state 1^(0)) to the state at times T, 2T, ■ ■ ■ . Equation (8) can be re- written as 

\z{T)) = exp (ii'yo-^ooT)] |^e-^("^+^)) . (10) 

Moreover, as the cyclic variation of a and P is continued, the state \n , a , P) picks up 
the same additional phase 7„ in each cycle - i.e., the same Hannay angle B is added in 
each interval T. Therefore, if 6k denotes the phase of the complex eigenvalue labelling 
\z{kT)) {k = 0, 1,2, ■■■), the evolution is equivalent to a rotation map [^] on a circle 
(S^), namely, 

4+1 = Ok- 2% A (mod 27r) (11) 

where 27rA = ujT + B. If A is an integer, every value of 9 is periodic, with period 
1. This is the case already discussed following Eq. (9). A more general possibility is 
A = p/q, a rational number. Again, every orbit of the map is periodic, but with a period 
q. Correspondingly, every initial state |-z(0)) has revivals at times gT, 2gT, ■ ■ ■ . But 
rational values of A constitute a set of measure zero. The generic case corresponds to 
an irrational value of A, and this is also the most interesting case. (It is also the most 
pertinent one from a practical point of view, as it relates directly to near revivals, as we 
shall see.) As is well known, the map no longer has any periodic orbits, but the iterates 
of any 9q cover S*^ densely as A; ^ 00. (Regarded as the Poincare section for motion on 
a 2-torus, this is the quasi-periodic case). The dynamics is ergodic but not mixing, with 



a uniform invariant density p{9) = l/(27r). Although we no longer have exact revivals in 
principle, the system comes arbitrarily close to these in practice (near revivals) - precisely 
the analog of Poincare recurrences. As a consequence, the following complete analysis of 
the statistics of near revivals becomes possible. 

Consider a prescribed small angular interval I^ of size 27re, located symetrically about 
the initial phase 6q on 5*^. Let A = [A] + 6, where [A] stands for the integer defined 
by A — 1 < [A] < A (for either sign of A), so that 6 is an irrational number satisfying 
< 6 < 1. It is then easy to see that 6^ E I^ ^ {^S} < e, where {x} denotes the fractional 
part of X. The corresponding correlation function then merely differs from unity by a term 
of order e^, because 

\{z{0)\z{kT))\'^ = exp [2|zp(cos27rA;A-l)" 

> exp [A\z\^TT^{k5y) > 1-A\z\^ttW (12) 

To order e, therefore, we may regard a return of 6^ to /^ as a (near) revival. The statistics 
of the occurrence times of such revivals is then identical to that of the recurrences to an 
angular interval of size 27re in the rotation map (11). The solution to the latter problem 



is given by certain gap theorems for interval exchange transformations [p2| , |23|| . Applying 
these to the case at hand, we obtain the following results in the long-time limit, after 
the transients due to specific initial conditions have died out and the invariant measure 
is attained. 

Ergodicity implies that the mean recurrence time (here, the mean time between suc- 
cessive near revivals) is, in units of T, the reciprocal of the invariant measure of I^. As this 
measure is uniform on the circle, the mean recurrence time is just T/e, as one may expect. 
The distribution of recurrence times is, however, quite remarkable |2^. In general (i.e.. 



for arbitrary e and 6), it is concentrated at no more than three points Ti = kiT, T2 = k2T 
and T3 = Ti + T2, where ki and ^2 are the least positive integers such that 

{ki5} < e and 1 - {M} < e (13) 

respectively. (Recall that e ^ 1; as long as e < 1/2, ki ^ k2.) It is easy to show that 
e < {ki6} + 1 — {k2S}: when the equality sign applies, only two recurrence times (Ti and 
T2) occur. If F{kT) denotes the normalized invariant (i.e., post-transient) probability 
that the revival time is kT {k = positive integer), then 



FikT) = (1/e) 



(e-{M})4,fci + (e-l + {M})'5, 



k,k2 



+ {{k,6} + l-{k26}-e)6kM^ 



k2 



(14) 



This is consistent with the requirement {k) = 1/e (which follows from ergodicity) provided 
/c2 {ki6} + ki{l — {k2S}) = 1, a relationship that can be established independently. With 
slight modifications, these results for near revivals continue to hold good in the periodic 
case {6 = rational number) as well. They are therefore generic, and do not require any 
fine-tuning of the parameters in the problem {u, T and B). 

Revivals are a manifestation of the interference arising from the different phases ac- 
quired under time evolution by the stationary basis states in the expansion of a non- 
stationary state. Hence it is the quantum-number-dependent part of the phase (both 
dynamical and geometric) that is relevant in this regard. That is why, in the present 
instance, it is only the coefficient B in the expression for 7„ (see Eq. (6)) that plays a 
role, and not the quantity A. (We have already seen why varying a does not lead to 
an n-dependent geometric phase, while varying f3 does so.) Therefore, as far as revivals 
are concerned, we may hold a fixed at the value zero throughout (i.e., displacement may 
be dispensed with altogether), without altering any of the conclusions. A cyclic varia- 
tion of the squeezing parameters jSi and (32 suffices to produce the relevant non-trivial 
anholonomy. For any given values of ujT (^ 1) and B, the near revivals of the corre- 
lation function concerned have been shown to be essentially equivalent to the Poincare 
recurrences in a rotation map on S*^. The steady-state distribution of recurrence times is 
explicitly determined. It is restricted to just three possible values. Their locations and 
relative frequencies of occurrence are easily found for any given 6 and prescribed e. 

As already stated, we have considered near revivals in the case of a linear spectrum. 
In the general nonlinear case (including, in the quantum optical context, multiphoton 



coherent states ||2^), additional interesting features appear, that have to be taken into 
account in establishing a mapping between near revivals (allowing for the effects of possible 
non-integrable phases) and Poincare recurrences. These include fractional revivals |^| and 



nonlinear dependence of the geometric phase on the quantum numbers. Details of this 
investigation will be reported elsewhere. 
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